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PROGRAM  SUMMARY 


This  report  summarizes  work  done  in  an  ongoing  program  under  contract 
RAD  80007  for  the  periods  October  1978  through  September  1979  (Project 
Order  RADC. 79-04)  and  October  1979  through  September  1980  (Project  Order 
RADC  80-01).  The  objective  is  to  gain  information  on  the  interaction  of 
ionizing  radiations  with  various  communications  materials.  Thus,  the 
intent  is  that  work  performed  under  this  contract  will  increase  our 
understanding  of  the  basic  mechanisms  involved  when  photons,  electrons, 
or  other  ionizing  radiations  interact  with  matter  in  its  condensed 
phases.  Emphasis  is  being  placed  on  studying  the  interactions  occurring 
in  those  materials  which  are  currently  used,  or  which  are  being  considered 
for  use,  in  electronics  systems.  Knowledge  of  quantities  such  as  the 
cross  sections  for  the  various  types  of  electron  interactions  in  these 
materials  enables  predictions  to  be  made  concerning  the  performance  of 
electronics  systems  in  various  situations,  such  as  exposure  to  the 
results  of  a  nuclear  explosion.  Our  approach  to  this  general  problem  is 
both  theoretical  and  experimental.  In  our  experimental  program,  we 
continue  to  make  observations  which  can  be  related  to  photon  and  electron 
interactions  with  solids.  The  theoretical  portions  of  the  program 
include  predictive  calculations  for  a  variety  of  interaction  phenomena 
as  well  as  interpretative  studies  in  support  of  our  experimental  efforts. 
Results  from  the  basic  theoretical  and  experimental  programs  are  employed 
in  applied  calculations  with  direct  relevance  to  the  performance  of 
electronics  systems  including  microsystems  produced  by  VLSI  technology. 

During  this  2-year  period  we  have  largely  concentrated  on  organic 
materials.  Optical  properties  have  been  measured  and  the  values  obtained 
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have  been  used  in  the  calculation  of  transport  parameters  for  ionizing 
radiations  incident  on  the  selected  organic  materials.  Our  studies  have 
led  us  to  propose  "universal"  relations  for  the  mean  free  paths  and 
stopping  powers  for  electrons  in  solid  organic  insulators.  We  have  also 
considered,  theoretically,  possible  mechanisms  associated  with  the 
protection  of  organic  crystals  against  radiation  damage.  Work  has 
continued  in  the  field  of  microdosimetry  with  the  calculation  of  both 
heavy-ion  track  structure  and  proton- induced  recoil  spectra  in  silicon. 

The  material  described  in  this  report  is  the  result  of  the  collab¬ 
orative  efforts  of  the  following  scientists: 

(1)  E.  T.  Arakawa,  J.  C.  Ashley,  R.  N.  Hamm,  R.  H.  Ritchie, 

J.  E.  Turner,  M.  W.  Williams,  and  H.  A.  Wright,  Health  and  Safety 
Research  Division,  and  V.  E.  Anderson,  Computer  Sciences  Division,  ORNL. 

(2)  C.  J.  Tung,  Institute  of  Nuclear  Science,  National  Tsing  Hua 
University,  Hsinchu,  Taiwan  300,  Republic  of  China. 

(3)  Martin  Pope,  Department  of  Chemistry,  New  York  University,  New 
York,  New  York  10003. 

(4)  T.  A.  Callcott  and  L.  R.  Painter,  Department  of  Physics  and 
Astronomy,  University  of  Tennessee,  Knoxville,  Tennessee  37916. 

(5)  M.  L.  Rustgi,  on  sabbatical  leave  from  the  Department  of 
Physics,  State  University  of  New  York  at  Buffalo,  Amherst,  New  York 
14260. 

Finally,  we  gratefully  acknowledge  the  continuing  advice  and  support  of 
J.  N.  Bradford  of  the  Rome  Air  Development  Center. 


I.  MEASUREMENT  OF  OPTICAL  PROPERTIES 


A  knowledge  of  the  optical  properties  of  a  material  over  an  energy 
range  which  covers  essentially  all  of  the  oscillator  strength  associated 
with  the  valence  electrons  enables  us  to  calculate  the  transport  parameters 
associated  with  the  interaction  of  ionizing  radiations  with  that  material. 
Thus,  for  a  number  of  years  we  have  been  measuring  the  optical  properties 
of  materials  over  an  extended  energy  range.  Over  the  past  few  years  we 
have  concentrated  on  organic  materials  in  a  condensed  phase.  Work  in 
this  reporting  period  has  involved  two  materials: 

A.  Polyethylene 

Work  started  in  the  previous  reporting  period  (see  report 
RADC-TR-79-234)  was  completed  and  published  [L.  R.  Painter,  E.  T.  Arakawa, 
M.  W.  Williams,  and  J.  C.  Ashley,  Radiat.  Res.  83,  1-18  (1980)]. 

B.  Kapton 

Measurements  were  initiated  and  completed  for  Kapton.  A  preprint 
of  the  results  from  this  work,  as  submitted  to  the  Journal  of  Applied 
Physios ,  forms  Appendix  A  of  this  report. 


II.  DETERMINATION  OF  ELECTRON  TRANSPORT  PARAMETERS 


A.  Mean  Excitation  Energies,  Mean  Free  Paths,  and  Stopping  Powers 

The  measured  optical  data  for  a  material  may  be  used  to  calculate 

various  electron  transport  parameters  for  that  material.  One  such 
parameter  is  the  mean  excitation  energy  I,  an  important  constant  in  the 
Bethe  theory  of  stopping  power.  Values  of  I  for  polystyrene  and  poly¬ 
ethylene  were  calculated  from  the  optical  data  for  these  solids.  We 
found  the  values  I  =  68.7  eV  for  polystyrene  [J.  C.  Ashley,  Phys.  Rev.  B 
19,  5429-5430  (1979)]  and  I  =  62  eV  for  polyethylene  [L.  R.  Painter, 

E.  T.  Arakawa,  M.  W.  Williams,  and  J.  C.  Ashley,  Radiat.  Res.  83,  1-18 
(1980)].  Other  parameters  which  can  be  calculated  from  a  knowledge  of 
the  optical  constants  over  an  extended  energy  range  are  the  mean  free 
path  and  the  energy  loss  per  unit  pathlength,  or  stopping  power,  for 
electrons  in  the  material.  A  model  insulator  theory,  with  parameters 
fixed  by  the  optical  data,  was  employed  to  calculate  inverse  mean  free 
paths  and  stopping  powers  for  low-energy  electrons  (10  eV  through  10  keV) 
in  polyethylene  [L.  R.  Painter  et  al.,  Radiat.  Res.  83,  1-18  (1980)]. 
Another  theoretical  model  involving  a  fit  to  the  energy-loss  function 
calculated  from  the  optical  data  was  used  to  estimate  electron  mean  free 
paths  and  stopping  powers  for  low-energy  electrons  in  Kapton  (see 
Appendix  A). 

B.  Ranges  of  Low-Energy  Electrons  in  Solids 

The  stopping  powers  of  Al ,  Au,  Ag,  and  Cu  for  electrons  were 
calculated  from  an  electron  gas  statistical  model  including  corrections 
to  the  stopping  power  due  to  the  if  effect.  From  these  stopping  powers 


the  continuous-slowing-down-approximation  (CSDA)  range  and  range  strag¬ 
gling  were  determined  for  electrons  with  energies  <10  keV.  Comparisons 
of  the  calculated  CSDA  ranges  with  experimental  data  were  made.  This 
work  was  published  in  the  following  reference:  C.  J.  Tung,  J.  C.  Ashley 
and  R.  H.  Ritchie,  IEEE  Trans,  mol.  Soi.  NS-26,  4874-4878  (1979). 
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III.  "UNIVERSAL"  CURVES 

A  concerted  effort  is  being  made  to  determine  "universal"  relations 
by  which  transport  parameters,  such  as  electron  mean  free  path,  can  be 
calculated  for  any  material  as  a  function  of  electron  energy.  During 
this  reporting  period  we  have  established  universal  relations  for  the 
mean  free  path  and  stopping  power  of  low-energy  electrons  (100  eV  to 
10  keV)  in  solid  organic  insulators. 

A.  Electron  Mean  Free  Paths 

A  simple  expression  involving  the  molecular  weight  and  density  of 
an  organic  sample  has  been  found  for  calculating  the  mean  free  path, 

A,  for  electrons  in  solid  organic  insulators  in  the  energy  range  100  eV 
to  10  keV.  The  expression  is  based  on  a  theoretical  formula  for  A  with 
input  parameters  derived  from  the  experimental  optical  spectra  for  solid 
organic  insulators  which  were  available  over  a  sufficiently  large  energy 
range.  These  data  consisted  of  the  optical  properties  of  DNA,  bovine 
plasma  albumin,  adenine,  guanine,  polyethylene,  and  polystyrene,  all  of 
which  had  been  obtained  in  our  laboratory  over  the  past  few  years.  As 
the  optical  properties  of  Kapton  (Appendix  A)  became  available,  it 
became  possible  to  check  our  universal  curve.  These  procedures  and 
results  are  outlined  in  the  following  submission  to  the  Proceedings  of 
the  Vlth  International  Conference  on  Vacuum  Ultraviolet  Radiation  Physics, 
University  of  Virginia,  Charlottesville,  Virginia,  June  2-6,  1980,  Vol. 
1-67,  pp.  1-3. 
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ELECTRON  MEAN  FREE  PATHS  IN  SOLID  ORGANIC  INSULATORS* 

M.  W.  Williams  and  J.  C.  Ashley 

Health  and  Safety  Research  Division 
Oak  Ridge  National  Laboratory 
Oak  Ridge,  Tennessee  37830 


The  concept  of  electron  mean  free  path,  a,  is  integral  to  calculations 
of  transport  phenomena  in  condensed  media.  In  the  case  of  biological 
molecules  this  parameter  is  difficult  to  measure  accurately,  so  we  have 
formulated  an  empirical  relation  from  which  the  mean  free  path  of  an  electron 
in  any  solid  organic  insulator  can  be  calculated  over  the  range  of  electron 
energies  from  100  eV  to  10  keV. 

The  mean  free  path  of  electrons  in  a  medium  can  be  calculated  from  the 
optical  properties  for  that  medium  over  the  whole  frequency  range.1  We 
calculated  A  as  a  function  of  electron  energy  L  from  100  eV  to  10  keV  for 
each  of  the  six  organic  insulators  (DNA,  bovine  plasma  albumin,  adenine, 
guanine,  polyethylene,  and  polystyrene)  for  which  the  optical  properties 
over  essentially  all  of  the  oscillator  strength  of  the  valence  electrons 
were  available.  For  each  material  the  relation  between  A  and  E  was  found 
to  be  of  the  form 


A(E)  =  E/[A  Lnl  ♦  B  -  C/E]  , 

where  the  constants  A,  B,  and  C  depend  on  the  material.  To  represent  the 
mean  free  paths  for  all  solid  organic  insulators  by  a  common  formula,  the 
assumption  was  then  made  that  A,  B,  and  C  are  each  functions  of  ny,  the 
valence  electron  density.  Analysis  then  showed  that  A»ny,  B«nv,  and  C«ny. 
Putting  our  results  in  terms  of  C  ■  pn/M,  where  v  is  the  density  of  the 
material  in  g^cm-’,  n  is  the  number  of  valence  electrons  per  molecule  or 
monomeric  unit,  and  M  is  the  molecular  weight  of  the  molecule  or  monomeric 


* 

Research  sponsored  jointly  by  the  Deputy  for  Electronic  Technology,  Air 
force  Systems  Command,  under  Interagency  Agreement  DOE  No.  40-226-70  and  the 
Office  of  Health  and  Environmental  Research,  U.S.  Department  of  Energy,  under 
contract  W-740B-eng-26  with  the  Union  Carbide  Corporation. 
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unit  in  g^mol'1,  we  obtained1  the  equation 

A£  =  E/[13.5  Ini  -  20.9  -  821/E]  ,  (1) 

which  was  predicted  to  relate  A  in  angstroms  to  £  in  electron  volts  for  any 
solid  organic  insulator  for  which  p,  n,  and  M  are  known. 

The  calculations  described  in  Ref,  1  have  been  repeated  with  a  restriction 
on  the  allowable  values  of  energy  transfer  to  account,  in  an  approximate 
fashion,  for  the  effects  of  exchange  between  the  incident  electron  and  the 
valence  electrons  in  the  solid.  In  addition,  an  increase  in  the  effective 
number  of  valence  electrons  by  a  factor  of  1.05  due  to  oscillator  strength 
coupling  between  core  and  valence  electrons  has  now  been  incorporated  in  the 
expression  for  With  these  changes,  but  still  neglecting  the  small  contri¬ 
butions  from  inner-shell  ionization,1  the  mean  free  path  may  now  be  written 

At  *  E/Cl 3.6  Lni  -  17.6  -  1400/E]  .  (2) 

A  plot  of  At  versus  E,  as  given  by  Eq.  2  is  shown  in  Fig.  1.  Mean  free  paths 
calculated  from  Eq.  2  are  '5%  smaller  than  those  predicted  by  Lq.  1  for  E.200  eV, 
equal  at  150  eV,  and  larger  at  100  eV.  Shown  on  either  side  of  this  curve 
are  the  plots  obtained  directly  from  the  optical  properties  for  the  individual 
materials,  bovine  plasma  albumin  and  polystyrene,  which  were  found  to  have 
the  largest  deviations  from  the  average  curve.  Based  on  these  data,  the 
universal  curve  given  by  Eq.  2  is  expected  to  predict  values  of  mean  free 
path  to  better  than  ♦ 1 0%  for  all  organic  insulators  over  the  energy  range 
from  100  eV  to  10  keV. 
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Fig.  1.  The  plot  of  At  versus 
E  as  given  by  Eq.  2  plus  the 
plots  for  bovine  plasma  albumin 
and  polystyrene  obtained  directly 
from  the  optical  properties. 


The  optical  properties  of  Kapton  Type  H  polyimide  film  have  been 
measured  recently  in  our  laboratory,  over  essentially  all  of  the  oscillator 
strength  of  the  valence  electrons.2  Figure  2  shows  the  values  of  A  from 
100  eV  to  10  keV  calculated  from  these  optical  data  and  as  predicted  from 
Eq.  2.  It  is  seen  that  there  is  excellent  agreement.  The  difference  between 
the  two  curves  ranges  from  5%  at  100  eV  to  02  at  10  keV.  This  is  well  within 
the  accuracy  claimed  for  the  universal  mean-free-path  relation  given  by  Eq.  2. 

It  is  anticipated  that  Eq.  2  should  also  predict  the  mean  free  paths 
of  electrons  in  liquid  organic  Insulators,  including  water.  Confirmation 
must  await  the  availability  of  optical  data  to  higher  energies. 


0UL*tMC  W-lOi* 


Fig.  2.  Comparison  of  the  mean 
free  path  of  Kapton  calculated  as 
a  function  of  electron  energy  from 
the  optical  data  and  from  the 
universal  relation  given  by  Eq.  2. 
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As  mentioned  at  the  end  of  the  previous  preprint,  we  anticipated 
that  this  approach  should  work  for  water.  Although  experimental  optical 
data  for  water  are  not  available  over  an  energy  range  that  covers  all 
the  oscillator  strength  of  the  valence  electrons,  an  estimated  spectrum 
exists  (R.  H.  Ritchie,  R.  N.  Hamm,  J.  E.  Turner,  and  H.  A.  Wright,  "The 
Interaction  of  Swift  Electrons  with  Liquid  Water,"  Vlth  Symposium  on 
Microdosimetry,  Eds.  J.  Booz  and  H.  6.  Ebert,  Brussels,  Belgium,  May  22- 
26,  1978,  Vol.  I,  pp.  345-354).  Shown  in  the  figure  below  is  the  dependence 
of  A  on  E  for  electrons  in  water  derived  from  the  universal  equation  for 
A  and  from  the  estimated  optical  spectrum  for  water.  It  is  seen  that 
there  is  excellent  agreement. 
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B.  Stopping  Powers  for  Electrons 

A  procedure  similar  to  that  described  in  part  A  has  been  used  to 
obtain  a  universal  relation  for  the  stopping  power  for  low-energy 
electrons  due  to  the  valence  electrons  in  organic  solids.  The  contri¬ 
butions  to  both  the  mean  free  path  and  the  stopping  power  due  to  K-shell 
ionization  in  carbon,  nitrogen,  and  oxygen  were  then  obtained  from 
theoretical  atomic  generalized  oscillator  strengths  and  added  to  the 
valence  electron  contributions.  This  work  is  described  in  the  preprint 
forming  Appendix  B.  This  paper  has  been  accepted  for  publication  in  the 
IEEE  Transactions  on  Nuclear  Science. 
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IV.  PROTECTION  AGAINST  RADIATION  DAMAGE  TO  ORGANIC  CRYSTALS: 

EXCITON  EFFECTS 

Experimental  measurements  by  Salih  and  Cosslett  [ Philos .  Mag.  30, 
225  (1974)]  indicate  that  radiation  damage  to  organic  crystals  of 
anthracene  and  coronene  may  be  reduced  by  factors  of  3  to  5  by  coating 

o 

the  crystals  with  thin  layers  (50-100  A)  of  aluminum  or  gold.  They 
suggested  that  these  results  might  be  explained  by  neutralization  of 
ions  or  holes  in  the  crystal  by  electrons  that  tunnel  from  the  metal 
overcoating. 

An  alternative  mechanism  to  explain  these  observations  has  been 
proposed,  i.e.,  that  excitons,  that  would  give  rise  to  lattice  defects 
in  the  absence  of  the  overcoating  layers,  are  quenched  at  the  metal - 
crystal  interface.  Quantitative  estimates  based  on  an  exciton  diffusion 
model  indicate  that  this  alternative  mechanism  may  be  important  in  the 
Salih-Cosslett  experiments.  Details  of,  and  predictions  from,  this 
model  have  been  published.  See:  R.  H.  Ritchie  and  Martin  Pope,  Philos. 
Mag.  B  40,  503-506  (1979). 


V.  MICRODOSIMETRY 


The  history  of  the  development  of  electronics  components  has  been 
one  of  continuous  reduction  in  size,  and  this  trend  is  expected  to  con¬ 
tinue  for  some  time.  In  present  day  large-scale  integrated  circuit 
technology,  critical  volumes,  such  as  the  active  region  in  a  memory 
chip,  have  been  reduced  to  dimensions  of  the  order  of  10  ym.  This  is  in 
the  size  domain  which  has  been  of  interest  in  the  field  of  microdosimetry 
in  connection  with  radiation  dose  deposition  of  significance  for  biological 
effects.  Microdosimetric  calculations  are  now  being  applied  to  help  in 
the  understanding  of  radiation  effects  in  solid-state  electronics  components 
In  the  past  few  years  the  phenomenon  of  cosmic-ray-induced  "soft 
errors"  in  satellite  computer  memories  has  been  observed.  This  is 
attributed  to  heavy  charged  particles  producing  sufficient  ionization 
within  a  critical  volume  to  change  the  bit  stored  in  that  location.  We 
have  undertaken  two  investigations  in  this  field.  In  one  we  have  used  a 
Monte  Carlo  transport  code  to  calculate  the  radial  dependence  of  dose 
deposited  in  the  vicinity  of  a  heavy-ion  track  in  silicon.  This  was 
undertaken  to  investigate  the  validity  of  error  induction  models  in 
which  it  is  assumed  that  a  track  may  be  treated  as  a  line  along  which 
energy  is  deposited.  In  the  other  we  have  calculated  the  spectrum  of 
heavy-recoil  particles  which  result  when  an  energetic  proton  collides 
with  a  silicon  nucleus. 

A.  Heavy-Ion  Track  Structure  in  Silicon 

The  energy  deposition  in  the  vicinity  of  a  heavy-ion  path  in  silicon 
was  investigated  by  a  Monte  Carlo  transport  analysis  of  the  delta  rays 
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produced  along  the  track.  The  dose  as  a  function  of  radial  distance  was 
calculated  for  C,  Al,  and  Fe  ions  with  energies  between  10  and  104  MeV. 

The  average  dose  deposited  in  a  cylinder,  parallel  to  the  track  and 
having  a  radius  between  10“3  and  1  pm,  as  a  function  of  the  distance  of 
the  axis  of  the  cylinder  from  the  path  was  also  determined. 

The  results  show  that  an  ion  track  is  very  small  in  lateral  extent 
as  compared  with  the  sizes  of  critical  regions  in  present  day  semiconductor 
memory  cells  (^10  pm).  This  remains  true  even  for  a  reduction  in  the 
dimensions  of  the  critical  region  by  an  order  of  magnitude.  It  is 
therefore  expected  that  a  model  in  which  ion  tracks  are  treated  as 
mathematical  lines  along  which  energy  is  deposited  in  accordance  with 
the  stopping  power  should  be  valid.  The  publication  describing  these 
calculations  and  results  is:  R.  N.  Hamm,  J.  E.  Turner,  H.  A.  Wright, 
and  R.  H.  Ritchie,  IEEE  Trans.  Hud.  Soi.  NS-26,  4892-4895  (1979). 

B.  Proton- Induced  Recoil  Spectra  in  Silicon 

Protons  ionize  too  sparsely  to  deposit  enough  charge  within  the 
critical  volumes  of  present  semiconductor  memory  chips  to  produce  soft 
errors.  However,  protons  may  undergo  nuclear  interactions  with  the 
nuclei  in  the  semiconductor  material  and  produce  heavy  recoils  which  are 
capable  of  depositing  a  sufficient  concentration  of  charge.  We  have 
used  a  nuclear  interaction  model  to  calculate  the  recoil  spectra  of 
heavy  fragments  produced  when  protons  are  incident  on  silicon. 

The  calculations  were  performed  using  MECC-7,  Medium-Energy  Intra¬ 
nuclear  Cascade  Code  System.  This  code  was  developed  at  Oak  Ridge 
National  Laboratory  and  is  available  from  the  Radiation  Shielding  Infor¬ 
mation  Center  at  Oak  Ridge. 
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In  this  code  an  incoming  nucleon  distributes  its  energy  among  its 
collision  partners,  the  nucleons,  inside  the  target  nucleus  until  it 
escapes  from  the  nucleus  or  until  it  loses  most  of  its  kinetic  energy 
and  at  last  is  captured  in  the  nuclear  potential.  The  evaporation 
process  which  follows  the  intranuclear  cascade  makes  use  of  the  work  of 
Le  Couteur  [Proc.  Phys.  Sos.  (London)  63A,  259  (1950)].  Use  is  made  of 
free  particle  nucleon-nucleon  data  whenever  cross-section  data  are 
needed.  These  data  as  well  as  the  nuclear  target  configurations  are 
specified  on  a  tape  used  by  MECC-7.  The  point  of  collision,  the  type  of 
collision,  the  momentum  of  the  struck  nucleon,  and  the  scattering  angles, 
etc.,  are  determined  by  statistical  sampling  techniques.  The  incident 
particles  may  be  protons  or  neutrons  with  energies  up  to  about  2500  MeV. 

The  option  exists  for  producing  at  least  one  pion  where  one  pion  production 
is  energetically  feasible. 

Calculations  were  made  for  protons  with  energies  of  20,  100,  130, 

200,  and  400  MeV  incident  on  28Si.  For  each  of  these  incident  proton 
energies,  we  obtained  the  cross  sections  for  production  of  various 
residual  nuclei  as  well  as  for  the  particles  directly  produced  by  the 
cascade  and  evaporation:  p,  n,  d,  t,  3He,  and  a.  For  each  of  these 
particles  and  residual  nuclei,  we  obtained  a  production  energy  spectrum. 

The  calculations  for  130  MeV  incident  protons  were  compared  with 
the  results  of  McNulty  et  al.  (P.  J.  McNulty,  G.  E.  Farrell,  R.  C.  Wyatt, 

P.  L.  Rothwell,  R.  C.  Filz,  and  J.  N.  Bradford,  to  be  published  in  IEEE 
Transactions  on  Nuclear  Science).  Though  the  fragment  mass  spectra  are 
very  similar,  two  differences  were  observed.  In  our  calculation  the 
mass  spectrum  peaked  at  the  mass  number  A=27,  while  in  McNulty  et  al.'s 
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calculation  it  peaked  at  A=25.  Secondly,  in  our  calculation,  no  frag¬ 
ments  with  5sAs9  were  observed.  In  the  following  figure  we  compare  our 
calculated  recoil  energy  spectra  with  their  measurement  of  energy  deposi¬ 
tion  in  a  silicon  slab  of  24.1-ym  thickness.  The  line  labeled  10<A<28 
is  our  calculated  absolute  cross  section  for  production  of  recoils  in 
the  indicated  mass  interval  having  energy  greater  than  E  as  a  function 
of  E.  The  straight-line  fit  to  their  experimental  data  is  normalized  to 
our  calculation  at  2  MeV.  The  other  line  is  our  calculated  results  when 
the  evaporated  particles,  p,  d,  t,  3He  and  a,  are  included.  It  is  also 
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normalized  to  the  other  curves  at  2  MeV.  We  find  an  exponential  dependence 
of  cross  section  on  recoil  energy  in  agreement  with  experiment,  but  our 
cross  section  falls  off  much  more  rapidly  with  increasing  energy,  and 
thus  we  would  predict  relatively  fewer  high-energy  recoils.  The  calculated 
cross  section  for  28Si  recoil  with  energy  greater  than  E  was  also  plotted 
but  could  not  be  compared  with  observations  for  lack  of  data. 

The  lack  of  agreement  could  be  attributed  to  the  exclusion  of 
several  important  processes  not  considered  in  the  code.  Among  the 
processes  which  become  important  at  higher  energies  is  the  spallation 
process  in  which  the  target  nucleus  is  spalled  into  several  pieces  by 
the  bombardment  of  energetic  protons.  Such  fragmentation  has  recently 
been  seen  when  target  nuclei  are  bombarded  with  a  particles,  12C,  160, 
and  other  heavier  projectiles.  The  energy  spectra  and  angular  distri¬ 
butions  of  these  fragments  differ  from  standard  compound-nucleus  theory 
expectations.  The  present  code  also  does  not  take  into  account  the 
knockout  of  certain  light  nuclei,  such  as  a  particles.  This  process  of 
course  becomes  important  at  higher  energies  while  the  compound  nucleus 
model  works  well  at  lower  energies. 
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ABSTRACT 

The  optical  properties  of  Kapton  Type  H  polyimide  film,  obtained  from 
reflectance  measurements,  are  presented  for  photon  energies  from  0.5  to  70  eV. 
Equations  for  the  energy-loss  function,  with  parameters  fixed  by  the  optical 
data,  are  employed  to  calculate  mean  free  paths  and  stopping  powers  for  low- 
energy  electrons  (100  eV  to  10  keV)  in  Kapton.  The  values  obtained  show  good 
agreement  with  those  predicted  by  recently  proposed  analytical  formulae  for 
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INTRODUCTION 


Direct  measurements  of  the  mean  free  paths  and  stopping  powers  of  low- 
energy  electrons  in  solid  organic  insulators,  such  as  polystyrene  and  poly¬ 
ethylene,  are  difficult.^  On  the  other  hand,  it  has  been  known  for  some  time 
that  a  knowledge  of  the  optical  properties  of  a  material  over  an  energy  range 
that  covers  essentially  all  of  the  oscillator  strength  of  the  valence  electrons, 

enables  the  transport  properties  for  ionizing  radiations  interacting  with  the 

2 

valence  electrons  to  be  calculated.  In  a  recent  publication  we  used  our  mea¬ 
sured  values  of  the  optical  properties  of  polyethylene  from  0.5  to  76  eV  in 
calculations  of  the  inelastic  mean  free  paths  and  stopping  powers  for  low-energy 
electrons  (10  eV  to  10  keV)  in  polyethylene.  Data  obtained  previously  on  five 
other  solid  organic  insulators  have  been  utilized  in  the  same  way^’^  and  simple, 
analytical  formulae  deduced  which  predict  mean  free  paths  and  stopping  powers 
for  electrons  in  solid  organic  insulators  in  the  energy  range  from  100  eV  to 
10  keV.  It  is  of  great  interest  to  check  these  predictions  using  a  material 
for  which  the  optical  properties  were  not  available  previously  over  an  energy 
range  which  covers  essentially  the  whole  of  the  oscillator  strength  for  the 
valence  electrons. 

Kapton  is  a  solid  organic  insulator  with  a  monomeric  unit  of  f  C22H10N2O5  4-. 

5 

Kapton  Type  H  polyimide  film  is  used  extensively  in  situations  demanding  good 
electrical  insulation,  radiation  resistance,  and/or  chemical  resistance,  particu¬ 
larly  at  high  operating  temperatures.  At  room  temperature  the  physical  properties 

5 

of  Kapton  are  similar  to  those  of  Mylar  polyester  film,  but  Kapton  is  less 


affected  by  temperature  changes  and  is  also  ^0  times  more  resistant  to 
radiation  damage  than  Mylar.  In  fact,  Kapton  has  been  found  to  retain  its 
useful  properties  over  the  wide  temperature  range  from  -269°C  to  +400°C.  In 


view  of  its  extensive  use,  the  optical  properties  of  Kapton,  and  hence  a  knowledge 
of  the  transport  parameters  for  ionizing  radiations  in  Kapton,  are  of  particular 
interest.  In  addition,  measurements  of  the  optical  properties  of  Kapton  over 
an  energy  range  which  covers  essentially  all  of  the  oscillator  strength  of  the 
valence  electrons  enables  a  comparison  to  be  made  between  the  values  of  electron 
mean  free  paths  and  stopping  powers  obtained  from  the  measured  optical  properties 
and  those  predicted  from  our  analytical  formulae.  Agreement,  within  the  limits 
of  experimental  errors,  would  lend  some  support  to  the  claim  that  our  analytical 
formulae  for  mean  free  path  and  stopping  power  apply  to  low-energy  electrons 
in  all  solid  organic  insulators. 

In  this  paper  we  present  the  optical  properties  of  Kapton  from  0.5  to  70  eV, 
and  then  use  these  data  to  calculate  electron  mean  free  paths  and  stopping 

powers  in  Kapton.  To  our  knowledge  the  only  optical  data  available  previously, 

5  6 

provided  by  the  manufacturer  and  appearing  in  the  open  literature,  were  in 

the  I.R.  and  in  the  vicinity  of  the  interband  absorption  edge. 

EXPERIMENTAL  PROCEDURE 

Kapton  is  available^  as  Type  H  film  with  nominal  thicknesses  from  7.6xl0-6  m 
to  1.27xl0"4  m.  Such  films  are  too  thick  for  transmission  measurements  in  the 
vacuum  UV  region.  Furthermore,  efforts  to  obtain  Kapton  films  suitable  for 
transmission  measurements  above  the  energy  region  of  the  UV  absorption  edge 
were  unsuccessful.  The  optical  properties  of  Kapton  were  thus  obtained  primarly 
from  reflectance  measurements,  using  photons  ranging  in  energy  from  0.5  eV 
(2500  nm)  to  70  eV  (17.7  nm)  reflected  from  Kapton  films  stretched  over  a 
sample  holder.  Measurements  were  made  with  a  spectrophotometer  (Cary  Model  14) 
in  the  0.5  to  5.4  eV  (2500  to  230  nm)  range,  a  Seya-Namioka  monochromator 
(McPherson  Model  235)  in  the  2.0  to  15  eV  (620  to  82.7  nm)  range,  and  with  a 


grazing  incidence  monochromator  (McPherson  Model  247)  in  the  13  to  70  eV 
(95.4  to  17.7  nm)  range.  The  values  obtained  for  the  normal  incidence 
reflectance  from  0.5  to  70  eV  are  shown  in  Fig.  1. 

For  photon  energies  where  the  extinction  coefficient,  k,  is  essentially 
zero,  interference  patterns  in  the  spectrum  obtained  by  reflection  were 
analyzed  to  give  the  refractive  index,  n,  and  the  thickness  of  the  film. 

Values  of  the  normal  incidence  reflectance  at  each  photon  energy  were  calcu¬ 
lated  from  n,  assuming  k  =  0.  For  photon  energies  where  k  is  small,  but  not 
zero,  reflectance  and  transmittance  measurements  were  analyzed  to  yield  n. 
Transmission  measurements  on  the  thinnest  Kapton  films  in  this  energy  region 
(i.e.,  in  the  region  of  the  ultraviolet  absorption  edge)  were  also  used  to 
obtain  the  values  of  k.  Values  of  the  normal  incidence  reflectance  were  then 
calculated  from  a  knowledge  of  n  and  k.  From  the  onset  of  absorption  to  ^25 
eV,  the  near-normal  incidence  reflectance  was  measured  directly  using  the  Seya- 
Namioka  and  grazing  incidence  monochromators.  In  the  3.0  to  5.4  eV  region, 
structure  in  the  normal  incidence  reflectance  was  confirmed  by  values  obtained 
from  reflection  measurements  in  the  Cary  14  spectrophotometer.  These  data 
were  normalized  at  3.0  eV.  For  photon  energies  *  25  eV,  at  each  energy  the 
measured  reflectance  as  a  function  of  the  angle  of  incidence  was  least-squares 
fit  to  Fresnel's  equations  to  yield  the  values  of  n  and  k.  For  photon  energies 
5  25  eV,  this  method  did  not  yield  accurate  values  of  n  and  k  for  Kapton.  In 
order  to  obtain  the  optical  properties  of  Kapton  over  the  whole  energy  range 
from  0.5  to  70  eV,  the  measured  n  and  k  values  from  0.5  to  3.0  eV  and  25  to  70 
eV  were  combined  with  the  values  of  the  near-normal  incidence  reflectance  from 
0.5  to  70  eV.  These  data  were  then  analyzed  by  a  Kramers-Kronig  technique 
which  we  had  previously  developed.^ 
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The  phase  change  on  reflection,  o,  is  given  as  a  function  of  incident 
photon  energy,  E,  by 
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7oeV 
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dE' 


dE' 
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where  -it  S  e  -  0.  The  integral  from  0  to  70  eV  was  calculated  as  a  function 
of  E  from  the  near-normal  incidence  reflectance  spectrum  shown  in  Fig.  1.  The 
integral  from  70  eV  to  infinite  energy  is  the  contribution,  ao(E),  to  the  phase 
change  on  reflection  due  to  absorption  in  the  energy  region  above  the  limit  of 
the  experimental  data.  Its  value  at  any  energy  E  depends  on  the  reflectance 
between  70  eV  and  infinite  energy.  For  energies  E  below  70  eV,  Ae( E )  must  be 
a  smoothly  varying  monotonic  function  of  E  with  d(Ae)/dEi0  and  d2(Ae)/dE2£0. 
Knowing  R  from  0  to  70  eV,  Ae(E)  can  be  calculated  from  Eq.  1  when  n  and  k  are 
known  since 


e(E)  =  tan-1[2k/(1-n2-k2)]  .  (2) 

The  calculated  values  of  ao,  from  0  to  3  eV  and  25  to  70  eV,  are  shown  in 
Fig.  2.  ao  was  interpolated  in  the  energy  range  from  3  to  25  eV  by  a  smoothly 
varying  function  as  shown  in  Fig.  2,  and  o  was  then  obtained  by  substitution 
back  into  Eq.  1.  Knowledge  of  o  and  R  in  the  interpolated  region  then  yielded 
n  and  k  from  the  equations 


n  =  (l-R)/(l+R-2\/R  coso) 


and  (3) 

k  =  (-2\/R  sino)/(l+R-2\/R  cose) 

This  method  of  Kramers-Kronig  analysis  does  not  require  extrapolation  of  the 
measured  reflectance  to  higher  energies.  The  values  of  n  and  k  obtained  from 
3  to  25  eV  are  normalized  to  directly  measured  values  in  the  adjoining  energy 
regions,  and  are  the  same  as  if  the  Kramers-Kronig  analysis  had  been  performed 
by  using  the  correct  extrapolation  of  R  to  infinity.  The  estimated  uncertainties 
on  the  values  obtained  for  n  and  k  are  everywhere  S  6%. 

DI SCU SSI  ON 

Our  results  for  Kapton  are  shown  in  Figs.  3  through  6,  in  the  form  of  n 
and  k,  c\  and  e2,  the  real  and  imaginary  parts  of  the  complex  dielectric 
function  e,  and  the  energy-loss  function,  Im(-l/e).  As  a  check  for  consistency 
of  these  data  with  the  results  for  other  hydrocarbons,  the  sum  rules  for  e2 
and  for  Im(-l/e)  are  shown  in  Fig.  7.  The  effective  number  of  electrons, 
neff,  per  monomeric  unit  of  Kapton  has  been  calculated  as  a  function  of  incident 

Q 

photon  energy  from  the  well-known  sum  rules, 

J  E'e2(E')dE'  =  — "eff 

o 

and 

f Vlm  (imf-'  •  — ¥-^  "eff  • 

o 


where  N0  is  the  number  of  monomeric  units  f  C22H10N205  )■  per  unit  volume  and 
e  and  m  are  the  charge  and  mass,  respectively,  of  an  electron.  By  comparison 
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2  9 

with  other  hydrocarbons,  ’  the  effective  number  of  electrons  participating 
up  to  70  eV  should  be  ^  93%  of  the  number  of  valence  electrons  in  the  mono¬ 
meric  unit.  For  Kapton  the  number  of  valence  electrons  is  138  and  hence  we 
anticipate  an  ng^  of  128  at  70  eV.  It  is  seen  that  ng^^  -  129  in  Fig.  7  -- 
well  within  the  experimental  uncertainty. 

APPLICATIONS 

A.  Electron  Inelastic  Mean  Free  Paths 

Values  of  inelastic  mean  free  paths  for  electrons  with  energies  from 

100  eV  to  10  keV  in  Kapton  were  calculated  from  the  optical  data  as  described 
3  10 

previously.  ’  Briefly,  a  model  for  the  energy-loss  function  was  employed 
which  requires  optical  data  to  fix  certain  adjustable  parameters.  The  values 
obtained  from  the  optical  data  are  shown  in  Fig.  8. 

We  have  predicted3, that  the  mean  free  path.  A,  for  solid  organic 
insulators  is  given  by 

1/A  =  (PN/M6Ml3.6Zn£-17.6-1400/rf)  ,  (4) 

O 

where  A  is  given  in  A  when  the  density  p  is  in  g/cm3,  the  molecular  weight  M 
is  in  g/mol,  the  electron  energy  £  is  in  eV,  and  N  is  the  number  of  valence 
electrons  per  molecule  or  monomeric  unit.  This  formula  was  obtained  by  calcu¬ 
lating  a  from  the  existing  optical  data  for  six  solid  organic  insulators  for 
which  the  optical  data  covered  essentially  the  whole  of  the  oscillator  strength 
of  the  valence  electrons.  Equation  4  was  then  the  best  fit  to  the  A  values  for 
the  six  materials.  This  has  been  termed  a  "universal"  curve  for  A  for  solid 
organic  insulators.  The  availability  of  the  optical  data  for  Kapton  over  a 
large  fraction  of  the  oscillator  strength  of  the  valence  electrons  gives  us 
an  opportunity  to  test  the  "universal"  curve  for  A.  Figure  8  shows  the  values 
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of  A  from  100  eV  to  10  keV  predicted  for  Kapton  from  Eq.  4.  It  is  seen  that 
there  is  excellent  agreement  between  the  two  curves  for  A.  The  difference 
ranges  from  5%  at  100  eV  to  <  3%  at  10  keV.  This  is  well  within  the  accuracy 
claimed  for  the  predictions  of  Eq.  4. 

It  is  anticipated  that  Eq.  4  should  also  predict  the  mean  free  paths  of 
electrons  in  liquid  organic  insulators  and  in  water.  Confirmation  must  await 
the  availability  of  optical  data  to  higher  energies. 

B.  Stopping  Powers 

In  a  similar  manner  the  stopping  power  for  low-energy  electrons  due  to 

interactions  with  the  valence  electrons  in  Kapton  can  be  calculated  from  the 

.  2  4 

measured  optical  properties.  ’ 

4 

It  has  been  predicted  that  the  contribution  of  the  valence  electrons, 

Sy,  to  the  stopping  power  of  a  solid  organic  insulator  is  given  by 

Sy[eV/A]  =  (pN/M£)  •(8.26xl02£«C -  2.61xl03  -  3.03xl04/fe)  ,  (5) 

for  100  eV  <  &  Z  10  keV,  in  terms  of  the  same  material  parameters  employed  in 
Eq.  (4). 

The  values  calculated  for  Sv,  both  from  the  optical  data  for  Kapton  and  from 
Eq.  (5)  are  shown  in  Fig.  9.  It  is  seen  that  there  is  satisfactory  agreement, 
with  9%  difference  at  100  eV,  3%  at  300  eV  and  less  than  1%  at  10  keV. 
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FIGURE  CAPTIONS 


Fig.  1  Normal  incidence  reflectance  of  Kapton  as  a  function  of  incident 
photon  energy. 

Fig.  2  Contribution  to  the  phase  angle  due  to  absorption  above  70  eV  as  a 
function  of  incident  photon  energy  for  Kapton  film. 

Fig.  3  The  optical  constants  n  and  k  of  Kapton  as  a  function  of  incident 
photon  energy. 

Fig.  4  The  extinction  coefficient,  k,  of  Kapton  as  a  function  of  incident 
photon  energy. 

Fig.  5  The  dielectric  functions,  ei  and  e2,  of  Kapton  as  a  function  of 
incident  photon  energy. 

Fig.  6  The  energy  loss  function  Im(-l/e)  of  Kapton  as  a  function  of  incident 
photon  energy. 

Fig.  7  The  effective  number  of  electrons  per  monomeric  unit  of  Kapton 
■f  C22H10N2O5  }  as  a  function  of  incident  photon  energy. 

Fig.  8  Comparison  of  the  mean  free  path  of  Kapton  calculated  as  a  function  of 
electron  energy  from  the  optical  data  and  from  the  universal  relation 
given  by  Eq.  4. 

Fig.  9  Comparison  of  the  valence  electron  contribution  to  the  stopping  power 
of  Kapton  calculated  as  a  function  of  electron  energy  from  the  optical 
data  and  from  the  universal  relation  given  by  Eq.  5. 
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ABSTRACT 

A  theoretical  description  of  the  inelastic  interactions  of  low- 
energy  electrons  (E^IO  keV)  with  solid  organic  insulators  is  presented. 
The  response  of  the  valence  electrons  to  energy  and  momentum  transfers 
is  inferred  from  a  model  for  the  energy-loss  function  with  parameters 
fixed  by  optical  data.  With  this  model  energy-loss  function,  inverse 
mean  free  path  and  stopping  power  due  to  interactions  with  the  valence 
electrons  are  calculated  for  several  organic  solids.  Contributions  to 
the  inverse  mean  free  path  and  stopping  power  due  to  K-shell  ionization 
in  carbon,  nitrogen,  and  oxygen  are  obtained  from  theoretical  atomic 
generalized  oscillator  strengths.  Simple,  analytical  formulae  are 
deduced  which  allow  predictions  of  mean  free  paths  and  stopping  powers 
for  electrons  in  the  energy  range  100  eV<E<10  keV  to  be  made  for  a 
large  group  of  organic  solids. 
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I.  INTRODUCTION 

Knowledge  of  the  inelastic  interactions  of  electrons  with  solids  is 
important  in  a  wide  variety  of  applications  which  involve  energy  deposition 
in  matter.  Information  generated  by  our  group2  on  electron  mean  free 
paths  and  stopping  powers  is  available  for  electron  energies  <10  keV 
for  several  metals  and  semiconductors  important  in  electronics  applications. 
Theoretical  calculations  of  these  transport  parameters  were  recently 
presented  for  the  organic  insulators  polystyrene3  and  polyethylene.4 
Such  material  by  material  studies  are,  however,  rather  time  consuming. 

It  would  be  useful  to  be  able  to  generalize  these  results  to  apply,  if 
only  approximately,  to  a  broad  group  of  materials.  The  purpose  of  this 
paper  is  to  describe  simple  analytical  formulae  for  mean  free  paths  and 
stopping  powers  for  electrons  of  energies  from  100  eV  to  10  keV  in  solid 
organic  insulators. 


II.  THEORETICAL  MODEL 

The  basic  equation  required  for  mean  free  path  or  stopping  power 
calculations  is  that  for  the  differential  inverse  mean  free  path  (DIMFP),^ 

infer  is?  r  ?  ,  o 

q- 

which  gives  the  probability  per  unit  distance  per  unit  energy  for  an 
electron  of  energy  E  to  suffer  an  energy  loss  in  the  interval  from  Hu  to 
Hu  +  d(Hu).  The  DIMFP  is  given  by  an  integral  over  possible  momentum 
transfers,  Hq,  to  a  medium  described  by  the  dielectric  response  function 
c(q,u).  The  medium  is  assumed  to  be  homogeneous  and  isotropic,  the 
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integration  limits  are  given  by  liq^  v7 2m[vT  +  v'E-'Fuo]  ,  and  aQ  h  Jf\2/me2 
is  the  Bohr  radius.  The  inverse  mean  free  path  p  and  stopping  power  S 
are  given  by 

0. 583  E 

m(E)  =  /  d(f)w)  dp/d(1iw)  , 

0 

and 

0. 583  E 

S(E)  =  /  ddto)  f)W  dji/d(-Piu>)  , 

0 

where  the  factor  0.583  in  the  upper  limit  of  these  integrals  accounts, 
in  an  approximate  fashion,  for  the  influence  of  exchange  between  the 
incident  electron  and  electrons  in  the  medium.6 

Since  the  dielectric  response  function  is  usually  not  known  in 
detail  for  a  given  material,  model  calculations  must  be  employed  to 
determine  e(q,w)  or,  as  required  in  Eq.  (1),  the  energy-loss  function 
Im[-l/e(q,w)].  The  model  used  here  involves  representing  Im[-l/e(q,w)] 
by  a  sum  of  Drude-type  functions  with  each  term  in  the  sum  involving 
three  adjustable  parameters.7’8  The  parameters  are  determined  for  a 
specific  material  by  fitting  the  model  energy-loss  function  for  q=0  to 
the  optical  energy  loss  function  lm[-l/e(0,w)]  as  determined  by  optical 
data  for  that  material.  Optical  data  for  photon  energies  up  to  ^80  eV 
for  polyethylene,  polystyrene,  DNA  (deoxyribonucleic  acid),  BPA  (bovine 
plasma  albumin),  adenine,  and  guanine  were  used  to  specify  the  model 
energy-loss  function  for  these  six  solid  organic  insulators.  Complete 
references  to  the  optical  data  are  given  in  Ref.  8.  The  energy-loss 
function  determined  in  this  manner  describes  only  the  response  of  the 
valence  electrons;  the  contributions  to  p  and  S  from  ionization  of  inner 
shells  is  treated  separately  and  will  be  discussed  below. 
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III.  MEAN  FREE  PATHS 

For  the  six  organic  insulators  for  which  sufficient  optical  data 
were  available,  inverse  mean  free  paths  were  calculated  from  Eqs.  (1) 
and  (2).  We  found8  that  the  contribution  to  the  inverse  mean  free  path 
due  to  interactions  with  the  valence  electrons  uv(E),  divided  by  the 
valence  electron  density  ny,  for  all  six  materials  could  be  represented 
to  better  than  ±10%  by  a  single  function  of  energy  for  100  eV<E^10  keV. 

The  calculations  described  in  Ref.  8  were  repeated  with  electron  exchange 
included  as  discussed  below  Eq.  (3).  In  addition,  an  increase  in  the 
effective  number  of  valence  electrons  by  a  factor  of  ^1.05  due  to 
oscillator  strength  coupling  between  core  and  valence  electrons,  as  noted 
earlier  in  polystyrene9  and  polyethylene,4  has  been  incorporated  in  the 
expressions  for  both  uv  and  Sy.  With  these  changes,  I  find 

Pv[A-l]/ny[au]  =  (152Z«E  -  198  -1 .57><104/E)/E  h  fy(E)  ,  (4) 

where  py  is  in  units  of  reciprocal  Angstroms,  E  is  in  electron  volts, 

and  ny  is  in  atomic  units.  For  simplicity,  the  factor  1.05,  discussed 
above  has  been  included  in  the  numerical  coefficients  in  Eq.  (4).  The 
use  of  this  factor  results  in  a  slight  overestimate  of  p  and  S  .  This 
point  is  discussed  further  in  Section  V.  v  v 

Values  of  ny  for  a  selection  of  organics  which  contain  only  nitrogen  and 
oxygen  in  addition  to  carbon  and  hydrogen  are  given  in  Table  I  along 
with  relevant  physical  properties  and  conversion  factors. 

Contributions  to  the  total  inverse  mean  free  path  due  to  ionization 
of  K  shells  in  carbon,  nitrogen,  and  oxygen  have  been  calculated  from 
theoretical  atomic  generalized  oscillator  strengths  (GOS's).  The  results 
for  carbon  and  oxygen  were  used  earlier  in  tables  of  electron  mean  free 


paths  in  polystyrene3  and  silicon  dioxide.10  The  values  for  nitrogen 
were  derived  from  McGuire's  GOS  values.11  The  ionization  cross  sections 
per  K  shell  o^,  in  atomic  units,  are  listed  in  Table  II  as  functions  of 
incident  electron  energy.  The  total  inverse  mean  free  path  is  given  as 
a  function  of  electron  energy  by  the  sum  of  the  valence  electron  contri¬ 
bution,  Eq.  (4),  and  K-shell  ionization  contribution  as 

m[A-  1 3  =  nv[au]fy(E)  +  1 .892oK^[au]nK^  .[au]  , 

where  the  sum  is  over  atomic  K  shells  for  carbon,  nitrogen,  and  oxygen 
and  n„  .  is  the  number  of  K  shells  per  unit  volume  (from  Table  I).  For 
each  material  in  Table  I  the  contribution  to  p  from  K-shell  ionization 
is  sl%.  Thus  electron  mean  free  paths  in  these  solids,  A=l/p,  are 
determined  mainly  by  interactions  with  the  valence  electrons,  Eq.  (4), 
and  the  inner-shell  contributions  may  be  omitted.  If  the  valence  electron 
density  is  expressed  in  terms  of  the  density  p  in  g/cm3,  the  number  of 
valence  electrons  per  molecule  or  monomeric  unit  N,  and  the  molecular 
weight  M  in  g/mol ,  the  inverse  mean  free  path  from  Eq.  (4)  may  be  written 

1/a[A]  =  (PN/ME).(13.6^E  -  17.6  -  1400/E)  , 

where  E  is  in  electron  volts.  This  equation  is  expected  to  predict 
values  of  electron  mean  free  paths  for  100  eV<E<10  keV  to  within  ±10% 
for  organic  insulators.  Mean  free  paths  calculated  from  Eq.  (6)  for 
E  s  200  eV  are  <3%  larger  than  those  predicted  earlier  for  electrons  in 
polystyrene, 3  and  <6%  larger  for  polyethylene,4  based  on  a  more  complicated 
model  for  the  energy-loss  function.  Comparison  of  our  earlier  theoretical 
result  with  the  meager  experimental  data  on  organics  is  given  in  Ref.  8. 
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Mean  free  paths  calculated  from  Eq.  (6)  are  <5%  smaller  than  those 
predicted  by  Eq.  (9)  of  Ref.  8  for  E  i200  eV. 

IV.  STOPPING  POWERS 

The  model  described  above  to  determine  the  energy-loss  function  has 
been  used  to  calculate  the  contribution  of  valence  electrons  to  the  mean 
energy  loss  per  unit  path  length  for  electrons  or  the  stopping  power  of 
the  solid  for  electrons.  I  find  that  contributions  to  the  stopping 
powers,  $y,  as  determined  by  calculations  for  polystyrene,  polyethylene, 
guanine,  and  adenine,  can  be  summarized  by  the  simple  result 

Sv[eV/A]  =  (PN/ME).(8.26xl027wE  -  2.61xl03  -  3.03xl04/E)  ,  (7) 

for  100  eVsE^lO  keV,  in  terms  of  the  same  material  parameters  employed 

in  Eq.  (6).  S  from  Eq.  (7)  may  slightly  overestimate  the  valence  electron 
contribution  to  stopping  power.  This  is  discussed  in  Section  V. 

Ionization  of  inner  shells  of  the  atoms  in  a  medium  can  make 
significant  contributions  to  the  stopping  power  for  electron  energies 
less  than  10  keV.  If  we  restrict  our  attention  to  organics  containing 
only  carbon,  nitrogen,  oxygen,  and  hydrogen,  a  group  which  includes  many  of 
the  important  materials  used  in  electronics  applications,  we  need  to 
consider  only  K-shell  ionization.  The  contributions  to  the  stopping 
power  due  to  K-shell  ionization  by  electrons  have  been  calculated  from 
GOS  values  as  described  previously.3  Values  of  sK,  the  stopping  power 
per  K  shell  per  unit  volume,  shown  in  Table  II,  were  extracted  from 
Ref.  3  for  carbon  and  from  Ref.  10  for  oxygen.  The  values  of  sK  for 
nitrogen  were  derived  from  McGuire's  GOS  tables.11 


Given  the  expression  for  Sy,  Eq.  (7),  and  . ,  where  the  subscript 
i  may  represent  carbon,  nitrogen,  or  oxygen,  the  total  stopping  power  is 
given  by 

S[eV/A]  =  Sv[eV/A]  +  51.4  -.[au^-Cau]  ,  (8) 

where  n^,  the  number  of  K  shells  per  unit  volume,  may  be  obtained  from 
Table  I.  The  numerical  factor  in  front  of  the  summation  sign  converts 

o 

the  sum,  expressed  in  atomic  units,  to  units  of  eV/A.  Stopping  powers 
calculated  from  Eq.  (8)  are  <5%  smaller  than  the  earlier  results  for 
polystyrene3  for  E  s  300  eV  and  ^15%  smaller  at  E  =  100  eV.  For  poly¬ 
ethylene4  Eq.  (8)  predicts  values  ^3%  smaller  at  E  =  10  keV,  increasing 
to  ^20%  smaller  at  100  eV,  than  our  previous  tabulated  values.  The 
differences  in  the  polyethylene  values  are  due,  in  part,  to  neglect  of 
exchange  in  our  earlier  tabulation.4  For  polyethylene  and  polystyrene, 

Eq.  (8)  predicts  stopping  powers  which  agree  within  2%  at  E  =  10  keV 
with  existing  tabulations12’13  for  E>10  keV. 

Equation  (8)  has  been  used  to  predict  stopping  powers  for  anthra¬ 
cene,  Bakelite,  Makrolon,  Plexiglas  (same  chemical  composition  as  Lucite, 
Perspex,  PMMA  resist).  Mylar,  cellophane,  and  two  types  of  nylon.  In 
all  cases  the  stopping  powers  at  10  keV  agree  with  tabulated  values12*13 
to  within  2%.  Stopping  powers  for  Kapton  have  also  been  predicted,  but 
no  tabulated  values  are  available  for  comparison.  The  K-shell  contributions 
to  the  stopping  powers  for  these  materials  accounted  for  slO-15%  of  the 
total  stopping  power  for  E  < 1 0  keV.  The  stopping  powers  of  some  of 
these  materials  are  displayed  in  Table  III.  It  is  interesting  that 
Kapton,  Mylar,  and  Bakelite,  which  are  quite  different  in  chemical 
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composition,  are  predicted  to  have  almost  identical  stopping  powers  for 
the  densities  shown  in  Table  I.  Note  also  that  S,  Eq.  (8),  is  linear  in 
p  and  thus  scales  easily  for  the  same  material  at  a  different  density. 

In  Fig.  1  we  show  S'hS/p  times  electron  energy  as  a  function  of 
electron  energy  for  Nylon  6,6,  Mylar,  and  Plexiglas.  The  solid  curves 
are  the  predictions  of  this  study;  the  dashed  curves  are  from  Brandt's 
tabulation. 14  Brandt's  results  agree  quite  well  with  other  tabulations1 2»1 3 
for  E>10  keV  and  fall  ^5-15%  below  my  predictions  at  E  =  256  eV. 

For  some  purposes  it  may  be  useful  to  have  a  simple  analytical 
expression  for  s^  rather  than  the  tabular  values  given  in  Table  II.  If 
E^  is  the  binding  energy  of  the  K-shell  electrons,  then  I  define  UK  by 
ILhE/E^.  The  values  for  s^[au]  in  Table  II  were  fit  by  an  expression  of 
the  form 

sK[au]EK[au]  =  (a  -  b/UK  +  c/U£)  .  (9) 

UK 

The  values  determined  for  a,  b,  c,  and  p  and  the  K-shell  binding  energies 
are  given  in  Table  IV.  As  shown  in  Table  IV,  Eq.  (9)  reproduces  the 
values  of  s^  from  Table  II  to  better  than  4%  for  s  1.5.  Since  the 
total  inner-shell  contribution  to  the  stopping  power  is  slO-15%  for  the 
materials  under  discussion  here,  the  use  of  Eq.  (9)  rather  than  sK  from 
Table  II  changes  the  total  stopping  power  by  <1%.  For  UK  <  1.5,  near 
threshold,  larger  deviations  from  the  values  in  Table  II  are  predicted 
by  Eq.  (9).  However,  the  values  of  sK  are  so  small  in  this  region  that 
even  a  factor  of  two  or  three  difference  produces  a  negligable  effect  on 
the  total  stopping  power. 
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It  is  known  that  cross  sections  for  K-shell  ionization  by  electrons 
satisfy  simple  scaling  laws14  (i.e.,  a^EjJ  should  be  a  universal  function 
of  U^).  Theory  indicates  that  a  similar  scaling  law  should  hold  for  s^ 
and  that  s^E^  should  be  a  function  only  of  for  all  atoms.  With  the 
values  of  s^  for  the  three  elements  from  Table  II  and  the  functional 
form  given  by  Eq.  (9),  I  find  that 

ln\S„ 

sK[au]EK[au]  =  (9.91  -  13.0/UK  +  3.13/U£)  ,  (10) 

K 

will  reproduce  the  Table  II  values  of  s^  for  carbon,  nitrogen,  and 
oxygen  to  within  7%  for  j  2.5,  within  20%  for  1.2  s  <  2.5,  and 
within  a  factor  of  ^2  for  Uj,  <  1.2.  The  discussion  following  Eq.  (9)  is 
alsr  appropriate  here;  the  use  of  Eq.  (10)  instead  of  the  values  from 
Table  II  should  lead  to  differences  of  <1%  in  the  total  stopping  power. 

V.  CONCLUSIONS 

I  have  presented  simple  analytical  expressions  for  predicting  mean 
free  paths  and  stopping  powers  for  electrons  in  organic  insulators  which 
should  be  useful  in  studies  of  energy  deposition  in  some  insulating 
materials  used  in  electronics  applications.  The  result  for  the  stopping 
power  may  not  be  restricted  to  insulators  since  the  presence  of  a  band 
gap  is  de-emphasized  by  the  term  tiw  in  the  integral  for  S  [see  Eq.  (3)]. 

In  addition,  all  the  formulae  may  apply  to  organic  liquids  as  well  as  solids. 

A  further  comment  is  nescessary  about  the  factor  1.05  included  in 
Eqs.  (4)  and  (7).  In  the  optical  limit,  q=0,  the  total  oscillator  strength 
for  K-shell  ionization  corresponds  to  less  than  two  electrons  per  K  shell 
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(1.72  for  C,  1.66  for  N  and  0).  The  valence  electron  density  was  increased 
by  a  factor  n  so  that  the  energy-loss  sum  rule,  du>  u>  Im[-l/e(q,a>)]  =  2*?n, 
yields  the  correct  total  electron  density  n.  For  q=0,  n:sl.05  for  the 
materials  in  Table  I.  If  n  were  independent  of  q,  then  pv  and  Sy  would 
be  increased  by  that  same  factor.  In  fact,  n  decreases  to  a  value  close 
to  one  for  large  q  so  a  different  correction  is  required  in  the  energy- 
loss  function  at  each  value  of  q.  The  results  for  m  and  S  would  then 

i  V  V 

be  multiplied  by  an  energy-dependent,  effective  value  of  n=neff(F)£l .05. 

Thus  the  stopping  powers  and  inverse  mean  free  paths  given  in  this  paper 
may  be  overestimated  by  a  few  percent,  with  larger  overestimates  occuring  at 
the  higher  energies.  This  point  and  a  more  detailed  discussion  of  stopping 
power  calculations  for  organic  insulators  and  the  extension  of  this  work 
to  organic  materials  other  than  solid  organic  insulators  will  be  presented 
elsewhere.16 
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Table  I.  Constants  for  Calculations  of  Mean  Free  Paths 
and  Stopping  Powers 


Organic 

Material 

b 

P 

Mc 

Nd 

n[au]xl04)e 

ny[au]xl02)^ 

c 

N 

0 

« 

- - - - 

Polyethylene 

4 

28.05 

12 

29.23 

3.507 

Polystyrene 

8 

1.05 

104.1 

40 

8.985 

3.594 

Adenine 

5 

5 

1.35 

135.1 

50 

8.904 

4.452 

Guanine 

5 

1 

5 

1.58 

151.1 

56 

9.317 

5.218 

Kapton 

22 

2 

5 

10 

1.42 

382.3 

138 

3.310 

4.568 

Anthracene 

14 

10 

1.25 

178.2 

66 

6.251 

4.125 

Bakelite 

44 

6 

36 

1.40 

660.7 

248 

1.888 

4.683 

Makrolon 

16 

3 

14 

1.20 

254.3 

96 

4.206 

4.038 

Plexiglas^ 

2 

8 

1.19 

100.1 

40 

10.59 

4.237 

Mylar*1 

4 

8 

1.385 

192.2 

72 

6.423 

4.625 

Cellophane 

5 

10 

1.50 

162.1 

64 

8.245 

5.277 

Rislan1 

22 

2 

2 

42 

1.045 

366.6 

152 

2.541 

3.862 

Nylon  6,6 

12 

2 

2 

22 

1.13 

226.3 

92 

4.450 

4.094 

Composition  of  material  in  terms  of  the  number  of  carbon  (C),  nitrogen 
(N),  oxygen  (0),  and  hydrogen  (H)  atoms  per  molecule  or  monomeric  unit. 

^Density  of  material  in  g/cm3. 

eMolecular  weight  in  g/mol  for  the  given  composition. 

^Number  of  valence  electrons  per  molecule  or  monomeric  unit. 

c>n[au]  is  the  number  of  molecules  or  monomeric  units  per  unit  volume 
in  atomic  units  (au);  n[au]  =  1 .480xl0-2!>  n[cm-3]. 

^nv[au]  is  the  valence  electron  density  in  atomic  units  (au); 
ny[au]  =  1.480xl0~25  ny[cnr3]. 

Colymethyl  methacrylate;  same  chemical  composition  as  Lucite, 

Perspex,  PMMA  resist. 

^Polyethylene  terephthalate;  same  chemical  composition  as  Melinex. 
tPolyamide  resin;  same  chemical  composition  as  Nylon  11. 


I 
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Table  II.  Contributions  to  Inverse  Mean  Free  Paths  and  Stopping 
Powers  Due  to  K-Shell  Ionization  in  Carbon  (C),  Nitrogen  (N), 

and  Oxygen  (0) 


o^taujxlO3 )a 

s 

K[au]xl02)Z; 

E  IeV] 

c 

N 

0 

c 

N 

0 

300 

0.131 

0.142 

400 

4.08 

4.93 

600 

9.66 

3.59 

0.163 

13.5 

5.95 

0.331 

800 

12.1 

5.69 

1.83 

18.9 

10.2 

4.03 

1000 

12.8 

6.70 

2.81 

22.0 

12.8 

6.62 

2000 

11.0 

6.77 

3.93 

23.1 

15.7 

11.3 

4000 

7.30 

5.03 

3.25 

17.2 

13.5 

11.1 

6000 

5.44 

3.92 

2.64 

12.8 

11.3 

9.54 

8000 

4.37 

3.29 

2.23 

10.3 

9.74 

8.32 

10000 

3.68 

2.78 

1.94 

8.67 

8.49 

7.39 

a 


a 


K[au]  is 


the  ionization  cross  section  per  K  shell  in  atomic  units. 


^sK[au]  is 
per  unit  volume 


the  stopping  power  due  to  K-shell  ionization  per  K 
in  atomic  units. 


shell 
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Table  Ill.  Stopping  Powers  (in  eV/A)  of  Some  Organic  Insulators 


E(eV) 

Kapton 

Mylar 

Bakel ite 

Plexiglas'3 

100 

4.57 

4.62 

4.67 

4.23 

150 

4.54 

4.59 

4.64 

4.20 

200 

4.14 

4.19 

4.24 

3.84 

300 

3.42 

3.46 

3.50 

3.17 

400 

2. 52 

2.95 

2.99 

2.70 

600 

2.30 

2.31 

2.36 

2.12 

800 

1.92 

1.93 

1.97 

1.76 

1,000 

1.66 

1.67 

1.71 

1.52 

2,000 

1.04 

1.04 

1.07 

0.944 

4,000 

0.621 

0.621 

0.636 

0.562 

6,000 

0.451 

0.450 

0.461 

0.407 

8,000 

0.357 

0.357 

0.365 

0.323 

10,000 

0.298 

0.298 

0. 304 

0.26S 

aPolymethyl  methacrylate;  also  Lucite,  Perspex,  PMMA  resist. 
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Table  IV.  Values  of  Parameters  Required  for  Predictions  of  Stopping  Power 
Per  K  Shell  Per  Unit  Volume  for  Electrons  with  E<10  keV  from  the  Equation 
s^CaujE^tau]  =  (a  -  b/U^  +  c/0^)Z.nUk/u£  for  Carbon  (C),  Nitrogen  (N),  and 
Oxygen  (0) 


c 

N 

0 

a 

16.4 

8.77 

10.1 

b 

38.1 

11.4 

15.9 

c 

26.6 

4.63 

7.98 

P 

1.15 

0.950 

1.00 

EK[auf 

10.4 

14.8 

19.7 

EK[eV] 

282 

402 

537 

%  Dlff.4 

$4% 

sU 

<2% 

aEK[au]  =  EK[eV]/27.2. 

This  is  the  percentage  difference  in  the  values 
predicted  by  the  fit  equation  for  U^sl.5,  using 
parameters  in  this  table,  from  those  presented 
in  Table  II. 


Fig.  1.  Stopping  powers,  S'=S/p,  as  functions  of  electron 
energy  for  three  solid  organic  insulators.  Solid  curves  from  this 
work;  dashed  curves  from  Ref.  14. 
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